Higher order coherence of exciton-polariton condensates 
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The second and third order coherence functions g^ n \0) (n — 2 and 3) of an exciton-polariton 
condensate is measured and compared to the theory. Contrary to an ideal photon laser, deviation 
from unity in the second and third order coherence functions is observed, thus showing a bunching 
effect, but not the characteristics of a standard thermal state with g^ n \0) = n\. The increase of 
bunching with the order of the coherence function, g^ 3 \0) > g^ 2 \0) > 1, indicates that the polariton 
condensate is different from coherent state, number state and thermal state. The experimental 
results are in agreement with the theoretical model where polariton-polariton and polariton-phonon 
interactions are responsible for the loss of temporal coherence. 

PACS numbers: 



In a seminal paper by R. J. Glauber marking the be- 
ginning of quantum optics, coherent states were studied 
as unique states, whose normalized coherence functions 
g( n )(0) are unity for all orders n It has been shown 
theoretically and experimentally that a photon laser is 
well characterized by coherence functions g^ n \0) close to 
one at pump intensities well above threshold, thus show- 
ing all order coherence, except for a random walk phase 
diffusion. Bose-Einstein condensation (BEC) of polari- 
tons has recently been the subject of detailed investiga- 
tions @, 0] and the question of its coherence properties 
arises. BEC is characterized by a large number of parti- 
cles sharing the same quantum state. Therefore it is ex- 
pected that the state of polariton condensates shows tem- 
poral coherence properties leading to a coherence func- 
tion of coherent state g( n \0) = 1 or TV-particle number 
state g( n \0) = 1 — However, there are important 

differences between photon laser and polariton BEC, one 
of which is that the former is based on stimulated emis- 
sion of photons from an excited level and population in- 
version is necessary to reach a threshold, while the lat- 
ter is based on stimulated cooling of polaritons to the 
ground state from excited states and BEC threshold con- 
dition is given by ~ 2 In Lj At, where n, At and L 
are the polariton density, thermal de Broglie wavelength 
and system size, respectively. Thus, when condensation 
occurs in a lower polariton (LP) ground state, the popu- 
lation inversion is not created yet in quantum wells 0, Q . 
Furthermore, scattering processes to the non-condensed 
polariton modes may influence the coherence properties 
of the polariton condensate. This last effect is absent 
in a photon laser. Contrary to a photon laser, devia- 
tion from unity in the second order coherence function 
g( 2 )(0) for the exciton-polariton condensate has been ob- 
served 0, 0] showing a bunching effect. The origin of 
the deviation from unity has not been well identified so 
far. However, in recent theoretical papers, the deviation 
from full coherence in higher order correlations have been 



shown [H, 0, • The model presented in @, 0, @] , consid- 
ers both polariton-polariton and polariton-phonon scat- 
tering. This scattering processes lead to two competing 
effects that involve the condensed polaritons: resonant 
scattering into and out of the polariton ground state 
involving one polariton in the ground state and para- 
metric scattering of two polaritons in the ground state 
into two states of opposite momentum. Below thresh- 
old, the statistics of the polaritons is determined by the 
thermal phonon bath, which is responsible for polariton 
relaxation, and the effect of polariton-polariton scatter- 
ing is negligible. In this case, the corresponding g( n \0) 
is given by that of a thermal state n\. Above thresh- 
old, resonant polariton-polariton and polariton-phonon 
interactions compete. The parametric scattering between 
ground state and excited states induces a depletion of 
the ground state population. This effect induces fluc- 
tuations in the population and manifests itself in a loss 
of coherence in the polariton condensate 0,0]. Indeed, 
when higher order coherence functions are calculated, a 
deviation from full coherence and bunching effect above 
threshold are obtained. 

Let us discuss this point in more detail. The particle 
statistics of a thermal state are characterized by the ge- 
ometrical distribution p(m) = j^(j^) m , where p(m) 
is a probability of m polaritons and \i is mean polari- 
ton number, and lead to a normalized coherence func- 
tion g^ n \0) = n!, which grows with n. In this case, we 
expect a larger bunching behavior in the higher order 
coherence functions, such as g( 3 \0) > g^ 2 \0). This be- 
havior is indeed found for polaritons in the ground state 
below threshold [Bj]. Above threshold, the calculated par- 
ticle distribution has a shape close to but different from 
a Poissonian distribution. The second order coherence 
function is larger than one and the higher order coherence 
functions are even larger [5]. In this paper, we present 
the measurement results of the second and third order 
normalized coherence functions. In particular, we show 
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that the observed behavior of growing bunching with the 
order n agrees well with the theory. 

The second order coherence function for exciton- 
polariton condensates has been measured for GaAs quan- 
tum well (QW) microcavity [2[ and CdTe QW microcav- 
ity @. Below condensation threshold, g^ 2 \0) takes al- 
most unity since both polariton emission decay time be- 
low threshold and detector response time are much longer 
than the intensity correlation time of polaritons. Above 
threshold g( 2 \0) increases steeply and then gradually 
saturates as the pump intensity goes up. A full coher- 
ence g^(0) = 1 has not been observed even at far above 
threshold. 

In the present work, we used 12 GaAs QWs embed- 
ded in GaAlAs/AlAs distributed Bragg reflector (DBR) 
microcavity that was need in our previous experiments 
@, [To[ . The sample has three stacks of four GaAs QWs 
which are embedded at the central three antinode po- 
sitions of a DBR planar microcavity. Lateral trapping 
potential is provided by a hole surrounded by a thin 
metal (Ti/Au) film which pushes photon- field amplitude 
antinode position at GaAlAs/air interface inside GaAlAs 
layer, resulting in the blue shift of the cavity resonance 
and hence LP energy under the metallic layer. By using 
this metallic hole structure, we could obtain the confine- 
ment of polaritons under a hole, which allows us to ac- 
cess to a single spatial mode in second and third order 
coherence measurements. An optical measurement sys- 
tem together with electronics for third order coherence 
measurement is shown in Fig. [TJ A mode-locked Ti:sa 
laser with 4 ps pulse width is utilized as a pump laser. It 
is focused on the sample surface with diameter of about 
50 fim. An objective lens collimates the output lumines- 
cence and near field image is made at the position of the 
pinhole by the second lens. Here we can spatially filter 
the metallic hole area of about 5 /am diameter. Then at 
the entrance slit of the spectrometer, far field is imaged 
by the third lens. Dispersion relation can be measured by 
the spectrometer and an attached nitrogen-cooled CCD 
camera 

The left-hand side of Figure[2]is an example of observed 
dispersion at 4P t h (Pth is threshold power ). In case we 
use the side exit of the monochromator, the output signal 
enters Hanbury-Brown Twiss (HBT) setup consisting of 
two non-polarizing beam splitters and photon detectors. 
We can measure (0) and (0) with this setup. @Sig- 
nal entering HBT setup is split into three paths by a 2:1 
splitting ratio beam splitter and a 1:1 ratio beam split- 
ter. We used single photon counting modules (Perkin 
Elmer, SPCM AQR series) for this coherence function 
measurement. The detectors' response time is relatively 
slow (about 300 ps timing resolution), but its slow re- 
sponse does not limit our capability of measuring g( 2 \0) 
and g( 3 \0) since we used a pulsed pump laser of 4 ps 
pulse width and the emission time becomes shorter than 
the correlation time at well above threshold. This setup 
allows us to pick up arbitrary part of condensates at LP 
ground state by setting appropriate wavelength filter of 




electronics 



FIG. 1: Experimental setup, pump: Mode- locked Ti:Sa laser 
with 76 MHz repetition rate and 4 ps pulse width. GaAs 
QW microcavity sample is hold in a cryostat which keeps the 
sample temperature at 10 K during experiment, ol is an ob- 
jective lens with 4 mm focal length and 0.55 NA. Momentum 
distribution is imaged at the entrance slit of the spectrome- 
ter by the ol and the following two lenses. Thus, dispersion 
curve of LP can be directly observed by the attached CCD 
camera and LP ground state emission is spectrally selected. 
When the signal goes through the other exit, it enters HBT 
setup. First two beam splitters (BS(2:1) and BS(1:1)) splits 
the signal into three paths with roughly equal intensities for 
maximizing 3-fold coincidence. Electronics consists of electri- 
cal processing units for taking 2-fold or 3-fold coincidences, i.e. 
single photon detectors (det 1 to det 3), delay units for timing 
adjustments, discriminator (delay units and discriminator are 
not shown in the figure), coincidence unit, and time- interval 
analyzer (TIA). 



spectrometer. The in-plane momentum region is deter- 
mined by collecting optics and actual CCD image size. 
Area surrounded by dotted line in Figure[2]is the detected 
region. 

For g^ (0) measurement, we need to take 3-fold coinci- 
dence of three detectors. For normalization of coherence 
function, 3-fold coincidence rate of simultaneous three 
photon detection for one pump pulse must be divided by 
accidental signal (product of independent signal rates of 
three detectors). In the actual measurement, we used a 
time interval analyzer (TIA) which measures time differ- 
ence between start and stop inputs. g^ 2 \0) can be di- 
rectly measured by TIA but g( 3 \0) cannot be since our 
TIA has just 2 input channels (start and stop). </ 3 )(0) is 
instead measured by the following method. At first, we 
took 2-fold coincidence between detector 1 and detector 2 
by 2-fold coincidence unit. Then the output signal is in- 
put to start port of the TIA (Fig. []]). Right-hand side of 
Figure [2] is an example of raw histogram taken by TIA. 
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FIG. 2: (left) Energy vs wavenumber dispersion of LP around 
4P t h- Closed area surrounded by dotted lines is picked up by 
using spectrometer and appropriate collection optics, (right) 
An example of 3 fold coincidence TIA data. Note that though 
it looks the same as usual (0) measurement, this histogram 
is triggered by input signals from 2- fold coincidence unit, 
therefore, this histogram gives information of # (3) (0)/# (2) (0). 
The highest peak corresponds to detections of three photons 
from one pump pulse. Surrounding peaks are 2-fold coinci- 
dence by one pump (start signal of TIA) but the stop signal 
of TIA is generated by other pump pulses at different time 
slots. Interval between each peak is 13 ns corresponding to 
the repetition rate of pump laser. 



The horizontal axis corresponds to delay time between 
start signal and stop signal. The interval between next 
pulses is 13 ns which corresponds to repetition period of 
the pump laser. The highest peak corresponds to 3-fold 
coincidence of simultaneously detected three photons for 
one pump pulse and surrounding peaks correspond to 2- 
fold coincidence by simultaneously detected two photons 
for one pump pulse and accidental third photon from dif- 
ferent time slot. So the surrounding peaks correspond to 
g( 2 )(0), i.e. just detector 1 and detector 2 are fired by 
the same pulse. In this configuration, we need to pre- 
liminary know g^ 2 \0) to obtain g^ 3 \0) since taking the 
ratio between central peak and average of surrounding 
peaks just gives g^ 3 \0) / g^ 2 \0) . So g^ 3 \0) can be ob- 
tained by multiplying g( 2 \0) to the ratio of central peak 
to surrounding peak. 

We measured g^ 2 \0) and g^ 3 \0) at various pump pow- 
ers (Fig. [3|). In Fig O theoretical results are drawn 
together with experimental data. The details of the 
theoretical model are given in Ref. [7]. Below thresh- 
old, though they are not shown in Fig 03 the statistics 
obeys thermal distribution, hence, the theory predicts 
g( n )(0) = n\. Just above threshold, they begin to gain 
coherence and rapidly decrease towards unity. However, 
the decreases stop around 1.5P t h clue to increasing ef- 
fect of polariton scattering. As the pump intensity in- 
creases, they converge into certain values (g^ 3 \0) ~ 2.5 
and g^(0) ~ 1.4). The experimental data shows g^ 2 \0) 
and g^ 3 \0) are still close to unity just above threshold. 
However, they increase and become closer to theoretical 
values as the pump intensity increases. After gradual 
increase, they reach a fiat area at around 4P t h- 
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FIG. 3: Pump intensity dependence of # (2) (0) and g { ' 6) (0). 
Pth is threshold power, square:*/ 2 ^), and circle:g3(0). Two 
curves are calculated coherence functions whose theoretical 
model is described in Q. 



To understand the discrepancies between theory and 
experiment at the pumping regime of 1 < P/Pth < 4, 
we need to consider intensity correlation time and decay 
time of the condensate emission. 

The ground state photoluminescence decay time mea- 
sured by a streak camera preceded by a spectrometer is 
shown in Fig. [H Since the streak camera was set after a 
spectrometer, we could avoid the contamination by the 
non-condensates and the wavelength range giving maxi- 
mum PL intensities at each pump power was picked up. 
Below threshold where photon statistics is expected to 
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FIG. 4: A photoluminescence decay time of a LP ground 
state. Around threshold, the decay time gradually decreases. 
It reaches a minimum value of 20 ps around 3-4 Pth where 
g^ 2 \0) and g^ 3 \0) approach to the theoretical values. 



obey thermal distribution, the observed statistics is al- 
ways g^(0) = g^(0) « 1 since decay time r of photolu- 
minescence is much longer than the intensity correlation 
time At c . In the case of longer r ^> Ar c , g^ n \r) = 1 
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0| . The measured is the integration of g^ n \r) over 
an integration time window which is given by an emis- 
sion time r in our case, , hence, it is close to unity after 
averaged over the whole emission lifetime. Just above 
threshold, the condition is still the same since the emis- 
sion pulse width is longer than the intensity correlation 
time until the pump rate reaches 4P t h- Finally far above 
threshold, the intensity correlation time becomes closer 
to the pulse width, and then the intrinsic noise property 
of the condensates begins to be detected properly. 

In conclusion, we have experimentally measured the 
second and third order coherence functions. The ob- 



served bunching effect is the experimental evidence for 
the relatively strong thermal and quantum depletion of 
the LP condensate. This higher order coherence function 
measurement technique may contribute to a further in- 
vestigation of coherence property of an exciton-polariton 
condensate. 
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